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We  describe  shape  transitions  of  thin,  solid  nematic 
sheets  with  smooth,  preprogrammed,  in-plane 
director  fields  patterned  across  the  surface  causing 
spatially  inhomogeneous  local  deformations.  A 
metric  description  of  the  local  deformations  is  used 
to  study  the  intrinsic  geometry  of  the  resulting 
surfaces  upon  exposure  to  stimuli  such  as  light  and 
heat.  We  highlight  specific  patterns  that  encode 
constant  Gaussian  curvature  of  prescribed  sign  and 
magnitude.  We  present  the  first  experimental  results 
for  such  programmed  solids,  and  they  qualitatively 
support  theory  for  both  positive  and  negative 
Gaussian  curvature  morphing  from  flat  sheets  on 
stimulation  by  light  or  heat.  We  review  logarithmic 
spiral  patterns  that  generate  cone /anti-cone  surfaces, 
and  introduce  spiral  director  fields  that  encode  non- 
localized  positive  and  negative  Gaussian  curvature 
on  punctured  discs,  including  spherical  caps  and 
spherical  spindles.  Conditions  are  derived  where 
these  cap-like,  photomechanically  responsive  regions 
can  be  anchored  in  inert  substrates  by  designing 
solutions  that  ensure  compatibility  with  the  geometric 
constraints  imposed  by  the  surrounding  media.  This 
integration  of  such  materials  is  a  precondition  for 
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their  exploitation  in  new  devices.  Finally,  we  consider  the  radial  extension  of  such  director 
fields  to  larger  sheets  using  nematic  textures  defined  on  annular  domains. 
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1.  Introduction 

It  is  well  known  that  inhomogeneous  local  deformations  in  thin  solid  sheets,  for  instance 
differential  growth  in  leaves,  can  lead  to  the  formation  of  Gaussian  curvature  and  complex  shape 
transitions  [1].  We  are  concerned  with  modern  responsive  materials  that  can  be  preprogrammed 
to  undergo  prescribed  spatially  inhomogeneous  expansions  and  contractions  in  response  to 
external  stimuli,  offering  exciting  possibilities  for  the  design  and  production  of  switchable 
surfaces  for  a  variety  of  applications  [2-6].  In  particular,  liquid  crystal  (LC)  solids,  either  glassy  or 
elastomeric,  have  orientational  order  about  a  director  n,  but  unlike  their  classical  fluid  analogues, 
they  cannot  flow.  Nematic  liquid  crystalline  glasses  and  elastomers  are  particularly  promising 
candidates  as  responsive  materials.  They  are  richer  than  isotropic  systems,  for  instance  those  with 
spatially  varying  capacity  to  swell,  in  that  both  the  direction  n  of  anisotropy  as  well  as  the  degree 
of  mechanical  response  are  variable.  Stimuli  inducing  such  changes  are  typically  light  or  heat, 
both  of  which  can  reduce  the  degree  of  order  and  cause  contraction.  The  contraction  X  <  1  along 
n,  and  elongation  by  X~v  >  1  in  the  directions  perpendicular  to  n,  are  large  in  elastomers,  where 
the  optothermal  Poisson  ratio  [7]  v  =  \  relates  the  perpendicular  and  parallel  responses  (with 
constancy  of  volume).  Glasses  have  X  ~  0.90-1.0  only,  and  their  v  2)  (volume  increase). 

Liquid  crystal  elastomers  have  a  subtle  mechanics,  because  director  rotation  can  accompany 
some  imposed  deformation  and  there  can  be  little  or  no  stress — a  remarkable  response  termed  soft 
elasticity.  Conceivably,  soft  elasticity  could  mask  the  Gaussian  curvature  we  explore,  except  that 
(i)  heating  all  the  way  to  the  isotropic  phase  eliminates  order  and  hence  anisotropy  that  can  be 
redirected  in  order  to  relieve  stress,  and  (ii)  we  find  states  that  are  unstretched  with  respect  to  the 
new  metric,  and  there  is  little  incentive  for  director  rotation.  We  henceforth  ignore  soft  elasticity 
and  treat  nematic  glasses  and  elastomers  in  the  same  way,  until  we  come  to  discuss  actuation  and 
device  implications. 

The  spontaneous  deformation  gradient  tensor  in  nematic  glasses  [8]  and  elastomers  [9]  is  of 
the  form 

F  =  (X  -  X~v)n  ®n  +  X~v  Id3,  (1.1) 

where  Id3  denotes  the  identity  operator  on  R3.  When  F  is  achieved,  the  heated  or  irradiated  body 
is  in  a  relaxed  state  without  stresses. 

The  distortions  associated  with  heating  a  sheet  with  uniform  n  are  quite  spectacular,  but 
more  subtle  and  perhaps  more  useful  effects  are  achieved  with  a  programmed-in,  non-uniform  n. 
Gaussian  curvature  arises,  because  the  metric,  specifying  the  intrinsic  geometry  of  the  deformed 
plane,  spatially  varies.  The  first  example  investigated  [4]  was  that  generated  by  concentric 
circles  of  n.  One  immediately  sees,  if  circumferences  contract  by  a  factor  of  X  <  1  and  radii 
extend  by  X~v  >  1,  the  circles  now  sit  on  what  has  morphed  into  a  cone,  concentric  with  its 
tip  where  the  localized  Gaussian  curvature  resides.  Experimentally,  effects  are  large,  even  in 
glasses  [10],  and  still  more  so  in  elastomers  [11].  These  director  circles  are  topological  charge  q  =  1 
defects  (disclinations)  in  a  two-dimensional  director  field.  Other  charges  instead  give  delocalized 
Gaussian  curvature  [7]  of  both  signs.  Programming  n  into  LC  sheets  is  possible  at  such  precision 
that  q  =  6  defects  can  be  created  in  elastomers,  and  in  the  case  of  glasses,  even  q  =  ±  10  defects 
have  been  produced  which  generate  highly  complex  surfaces  [11,12].  In  optics,  such  structures  are 
used  as  axially  symmetric  wave  plates  [13],  for  instance  with  charge  q  =  64,  or  termed  '^-plates', 
producing  sophisticated  light  polarization,  with  charge  q  =  100  being  employed  [14]. 

More  general  distributions  of  Gaussian  curvature  by  variation  of  director  in  plane  were 
attacked  in  the  seminal  work  of  Aharoni  et  al.  [5]  that  describes  the  interplay  between  the  nematic 
director  field  of  a  thin  elastomeric  sheet  and  the  resulting  three-dimensional  configuration 

attained  upon  heating.  In  particular,  they  consider  the  difficult  reverse  problem  of  constructing 
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a  director  field  that  induces  a  specified  two-dimensional  metric.  In  cases  of  one-dimensional 
variation  (leading  to  surfaces  of  revolution),  they  have  a  direct  way  of  finding  the  metric.  These 
authors  also  go  on  to  consider  more  ambitious  problems  such  as  two-dimensional  variation. 
In  this  paper,  we  follow  [6]  and  consider  some  specific  in-plane  director  field  patterns  on  thin 
nematic  sheets  that  give  interesting  curvature  distributions  and  consequent  shapes  that  are 
then  compared  qualitatively  with  our  experiments.  We  continue  by  considering  spiral  director 
patterns  that  produce  more  advanced  shapes — cones,  hyperbolic  cones,  (pseudo)spherical  caps 
and  spindles.  An  additional  but  important  constraint  is  then  added  in — that  the  generated  shape 
is  geometrically  compatible  with  a  surrounding,  inert  sheet,  so  that  the  morphing  structure 
can  be  anchored.  Practical  devices  will  require  such  attachment,  as  they  respond  to  stimuli,  in 
order  to  exploit  induced  Gaussian  curvature  to  pump,  do  work  or  actuate.  We  will  conclude  by 
discussing  how  applications  might  be  achieved  because  of  this  essential  advance.  We  now  review 
the  mathematical  framework  [5,6]  used  in  the  subsequent  analysis. 
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2.  Generating  intrinsically  curved  surfaces 

(a)  Gaussian  curvature  through  patterning 

It  is  assumed  that  the  director  field  does  not  vary  across  the  thickness  of  the  sheet,  so  that  the 
same  pattern  is  repeated  at  each  level  of  thickness.  For  sufficiently  thin  sheets,  stimulation  of 
the  system  will  result  in  pure  bending  of  the  sheet  at  no  stretch  energy  cost  and  one  expects 
an  isometric  immersion  of  the  prescribed  local  deformations  as  determined  by  the  director  field 
pattern. 

Let  (xi,X2)  £  co  c  R2 3  be  Cartesian  coordinates  parametrizing  the  midsurface  of  the  initially  flat 
sheet  and  n(x  1,^2)  =  n\el  +  be  the  director  field  pattern  across  the  surface,  where  el ,  e 2  form 
the  standard  orthonormal  basis  of  R2.  The  associated  in-plane  spontaneous  deformation  tensor  F 
has  components  Fap  =  (k  —  k~v)nanp  +  k~v8ap,  where  a,  =  1,2.  The  resulting  two-dimensional 
metric  of  the  deformed  sheet  upon  stimulation  is  a  =  FTF,  which  simplifies  to 

aap  =  (X2 -X~2v)nanp+k~2vSap.  (2.1) 

The  Gaussian  curvature  K  of  a  surface  is  an  intrinsic  geometric  property  ( Theorema  Egregium)  that 
is  determined  by  the  first  fundamental  form  aap  of  the  surface  via 

K  =  -T^r2 - a2r2  +  r^r2  - r^r2  +  r2r2  - r2r2),  (2.2) 

an 

where  F^p  =  j ttT(r  (dpL^op  T  dpftpia  ~  and  ( a ^°)  =  (#pp)  ^ • 

Ignoring  spatial  variation  in  A  that  would  also  contribute  to  the  Christoffel  symbols  r*p  and 
hence  K}  we  concentrate  on  the  two-dimensional  director  field  that  is  characterized  by  an  angle 
scalar  field  x/r  =  ^(^1,^2)  specifying  the  in-plane  orientation  of  the  director  at  each  point  on  the 
initially  flat  sheet,  so  that  ni  =  cos  x/r  and  ft  2  =  sin  x/f.  The  Gaussian  curvature  determined  by  the 
nematic  metric  can  be  expressed  in  terms  of  x/r  as 

x2v  _  i-2 

K  = - - - [(dfa  ~  ^1  ^  —  431^32^)  sin(2^)  +  2(3i32^  +  (32 VO2  —  (9iV02)  cos(2^)].  (2.3) 

We  note  here  that  if  we  rotate  the  director  associated  with  a  given  pattern  by  7t/2  in-plane  at  every 
point,  so  that  x/r  ^  +  tt/2,  then  the  resulting  Gaussian  curvature  flips  sign  everywhere,  because 


1  Indeed,  k(r)  variation,  but  in  isotropic  systems  without  a  director,  is  the  basis  of  Gaussian  curvature  by  spatially  dependent 

swelling  [2,3]. 

3 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


rspa.royalsocietypublishing.org  Proc.  R.  Soc.  A  472: 20160112 


Downloaded  from  http://rspa.royalsocietypublishing.org/on  March  10,  2017 


sin  2^  — ►  —  sin  2^  and  cos  2^  — ►  —  cos  2^.  That  is, 

K^-K  (2-4) 

We  refer  to  a  pair  of  director  field  patterns  that  are  related  by  a  7i/2  radian  rotation  of  the  directors 
as  orthogonal  duals,  pairs  that  relate  to  heating  and  cooling.2  In  this  paper  we  consider  only  X  <  1, 
that  is  heating;  we  achieve  changes  in  sign  of  K  by  changing  the  n(r). 

(b)  Spherical  and  pseudo-spherical  surfaces 

We  now  restrict  our  attention  to  director  fields  of  the  form  tt  =  cos  i/r(x 2)  e1  +  sin  ^(*2)  c2,  whose 
alignment  angle  field  varies  only  with  respect  to  one  of  the  coordinates.  The  Gaussian  curvature 
upon  stimulation  is  K  =  —  \{X~2  —  A,2v)(y V"  smlx/f  +  2\j/'2  cos  2i/r).  We  can  rewrite  this  as 

d2 

-j— 2  cos  2^r  =  4  C(K),  (2.5) 

where  C(K)  =  K/(X~2  —  X2v),  and  solve  for  constant  K  >  0  to  find 

t/c(x2)  =  ±2  arccos(ci  +  C2X2  +  2C(K)X2)/  (2.6) 

where  c\,  c2  are  constants  of  integration.  This  pattern  generates  constant  Gaussian  curvature  K 
wherever  it  is  well  defined.  Now,  consider  the  particular  solution 

^(*2)  =  arccos  ^2  (l  —  ^r)  -  lj  ,  (2.7) 

corresponding  to  c\  =  1,  c2  =  — 4/L,  and 

1  ■ -2  _  i2v 

— -  >0'  (2-8) 

This  can  be  rewritten  as  t/c(x2)  =  arccos (1  —  x2/L),  which  describes  a  well-defined  pattern  for 
0  <  x2  <  2 L.  The  pattern  on  the  square  domain  =  [0, 2L]  x  [0, 2L]  is  shown  in  figure  la.  By 
integrating  along  the  director  field  lines,  we  note  that  the  integral  curves  of  this  pattern  consist  of 
semicircles  of  radius  L  that  are  shifted  along  the  xi-axis.  Thus,  if  we  seek  to  encode  a  particular 
constant  positive  Gaussian  curvature  K  =  Kq  >  0  across  an  initially  flat  sheet,  then  we  can  do  so  by 
encoding  the  pattern  obtained  by  shifting  a  semicircle  of  radius  L  =  (1  /s/Kq)^/x~2  —  X2v  as  shown 
in  figure  lb. 

By  the  observation  that  the  orthogonal  dual  of  a  given  two-dimensional  director  field  pattern 
generates  the  exact  opposite  Gaussian  curvature  at  every  point,  we  can  encode  constant  negative 
Gaussian  curvature  K  =  —Kq  on  a  thin  nematic  sheet  by  simply  using  the  orthogonal  dual 
of  the  pattern  that  encodes  positive  curvature  Kq  >  0.  Returning  to  the  example  of  figure  la, 
where  a  pattern  encoding  constant  positive  curvature  K  =  (X~2  —  X2v)/L2  was  defined  on  the 
square  domain  co  =  [0,2L]  x  [0,2L]  by  ^(x2)  =  arccos(l  —  x2/L),  we  immediately  obtain  a  pattern 
on  the  same  domain  which  encodes  constant  negative  Gaussian  curvature  K  =  —  (X~2  —  X2v)/L2, 
by  simply  taking  \jr(x 2)  =  arccos(l  —  x2/L)  +  7i/2.  The  resulting  pattern  is  shown  in  figure  lc.  This 
pattern  is  generated  by  shifting  a  tractrix  curve  along  its  axis. 

For  a  surface  in  R3,  the  components  aap  and  bap  of  the  first  and  second  fundamental  forms 
satisfy  a  system  of  algebraic  differential  equations  known  as  the  Gauss-Codazzi-Mainardi 
equations.  Conversely,  any  pair  (a,  b)  consisting  of  a  symmetric  and  positive  definite  matrix  field 
(aap)  and  a  symmetric  matrix  field  {bap)  that  satisfy  the  Gauss-Codazzi-Mainardi  equations 
determines  a  unique  surface  up  to  a  rigid  transformation  in  R3  [15].  Thus,  to  determine  the 
equilibrium  configuration  of  the  mid-surface  of  an  initially  flat  nematic  sheet  upon  stimulation, 
we  also  need  to  know  the  components  bap  of  the  second  fundamental  form  that  minimize  the 
bending  energy  subject  to  the  Gauss-Codazzi-Mainardi  constraints. 

2 Heating  the  orthogonal  dual  to  give  a  k  <  1  along  its  director  corresponds  to  cooling  the  original  pattern  to  give  an  elongation 
along  its  director  of  A  =  X~v  >  1.  Thus,  heating  and  cooling  give  converse  Gaussian  curvatures. 
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Figure  1.  (a)  The  director  field  defined  by  i/r(x2)  =  arccos(l  —  x2/L)  on  the  square  domain  co  =  [0,2 L]  x 
[0,2L].  This  pattern  generates  constant  Gaussian  curvature  K  —  (A-2  —  A2y)/L2  >  0  upon  stimulation,  (b)  The 
nematic  pattern  obtained  by  shifting  a  semicircular  arc  of  radius  L  =  (1  /\/K)Va-2  —  A2y  along  the  x\  -axis  generates 
constant  positive  Gaussian  curvature  K  >  0  upon  stimulation,  (c)  The  director  field  defined  by  if  (x2)  =  n /2  +  arccos 
(1  —  X2/L)  on  the  square  domain  co  =  [0, 2 L]  x  [0, 2 L].  This  pattern  generates  constant  negative  Gaussian  curvature 
K=  —(A-2  —  A 2y)/L2  <  0.  (Online  version  in  colour.) 


For  a  fixed  two-dimensional  metric,  the  problem  of  identifying  equilibrium  configurations  that 
minimize  the  bending  energy  reduces  to  the  problem  of  minimizing  the  Willmore  functional 


Iw  = 


H2  dS, 


(2.9) 


where  H  is  the  mean  curvature  of  the  deformed  surface,  among  isometric  immersions  of  the  given 
metric  [16-18]. 

For  a  metric  of  constant  positive  Gaussian  curvature,  it  is  easy  to  show  that  the  Willmore 
functional  is  minimized  precisely  for  spherical  solutions.  That  is,  a  flat  nematic  sheet  whose 
director  field  encodes  constant  positive  curvature  K  is  expected  to  form  part  of  a  sphere  of  radius 
R  =  1/VK  upon  stimulation,  assuming  that  the  sheet  is  small  enough  to  exclude  the  possibility  of 
self-intersection  [6]. 

In  the  case  of  a  metric  of  constant  negative  Gaussian  curvature,  identifying  minimizers  of 
the  Willmore  functional  is  considerably  less  straightforward.  In  [19],  it  is  shown  that  for  a 
hyperbolic  elastic  disc  that  has  already  undergone  local  deformations,  surfaces  that  are  geodesic 
discs  lying  on  hyperboloids  of  revolution  of  constant  Gaussian  curvature  are  minimizers  of  the 
Willmore  functional  among  smooth  immersions  of  the  metric.  These  solutions  will  appear  as 
saddle  shapes  in  experiments  and  are  expected  to  be  energetically  favourable  for  sufficiently  small 
discs.  However,  it  has  been  shown  numerically  that  certain  non-smooth  wavy  surfaces  formed  as 
odd  periodic  extensions  of  subsets  of  so-called  Amsler  surfaces  are  energetically  more  favourable 
than  the  smooth  saddle  shapes  that  correspond  to  discs  lying  on  hyperboloids  of  revolution  when 
the  radius  of  the  hyperbolic  disc  is  sufficiently  large  [20].  We  synthesize  LC  solid  films  with 
spatially  programmed  directors  in  order  to  realize  shape-changing  surfaces  of  constant  Gaussian 
curvature. 


(c)  Experimental  investigations 

The  director  profile  in  the  plane  of  nematic  LC  solid  sheets  can  be  programmed  through  a  variety 
of  methods,  including  mechanical  and  magnetic  fields  [21,22].  Using  these  methods,  however, 
it  is  difficult  to  spatially  control  the  director  orientation.  Here  we  use  chemistry  amenable  to 
surface  alignment:  precursor  molecules,  to  what  will  become  the  nematic  solid  sheet,  have  been 
specifically  designed  to  align  to  treated  surfaces.  Using  this  approach,  low  molar  mass  nematic  LC 
monomers  are  filled  between  two  plates  separated  by  a  well-defined  gap.  The  treated  surfaces, 

on  the  interior  of  the  plates,  direct  the  ordering  of  the  LCs  along  a  specific  orientation  through 
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Figure  2.  (a)  Polarized  optical  images  of  the  patterned  director  profiles  predicted  to  generate  positive  (left)  and  negative  (right) 
curvature.  The  patterns  are  optically  equivalent  between  crossed  polarizers.  The  director  orientation  at  the  edges  and  centre  of 
the  pattern  is  indicated  with  arrows.  Each  square  film  has  a  side  length  of  10  mm.  [b)  Positive  (left)  and  negative  (right)  Gaussian 
curvature  in  15  jam  thick  glassy  LC  solid  film  at  175°C.  (Online  version  in  colour.) 


Figure  3.  From  left  to  right:  (a)  the  initially  flat  configuration  of  a  circular  glassy  film  15  jam  in  thickness  and  7.1  mm  in  diameter. 
(b)  The  positive  Gaussian  curvature  pattern.  The  dashed  circle  indicates  the  boundaries  of  the  circular  film,  (c)  The  formation  of 
positive  Gaussian  curvature  in  the  actuated  state  from  two  distinct  viewing  angles,  (d)  The  negative  curvature  pattern  obtained 
as  the  orthogonal  dual  director  field,  (e)  The  formation  of  negative  Gaussian  curvature  in  the  actuated  state  from  two  viewing 
angles.  (Online  version  in  colour.) 


the  thickness  of  the  material.  By  using  reactive  nematic  mesogens,  this  director  orientation  can 
be  trapped  in  an  elastic  solid.  Spatially  complex  director  patterns  in  LC  cells  are  prepared  using 
point-by-point  photoalignment  of  an  azobenzene  dye  by  irradiation  with  polarized  light  [12,23]. 
By  altering  the  polarization  of  the  incident  light,  the  in-plane  orientation  of  the  director  of  the  LC 
can  be  spatially  controlled.  The  resulting  director  field  is  a  pixelated  approximation  of  the  desired 
smooth  pattern  with  each  pixel  measuring  100  x  100  pm. 

A  number  of  glassy  liquid  crystalline  solids  have  been  demonstrated  to  be  compatible  with 
surface  alignment  techniques.  Here  we  use  one  such  composition  with  X  =  0.94  and  v  =  0.92 
[24].  Specifically,  we  use  the  composition  with  the  lowest  cross-link  density  from  this  work.  This 
composition  is  representative  of  the  larger  class  of  nematic  LC  glasses  that  can  be  aligned  using 
surface  alignment  techniques  [25].  The  director  patterns  depicted  in  figure  1  a,c  were  chosen  to 
assess  the  viability  of  generating  Gaussian  curvature  on  exposure  to  stimulus.  After  fabrication, 
the  LC  solid  film  is  flat  at  25° C  and  retains  the  expected  birefringence  of  an  aligned  nematic,  as 
seen  in  figure  2a. 

The  thermally  induced  shape  change  of  the  nematic  LC  glass  is  shown  in  figure  2b.  As 
predicted,  the  pattern  depicted  in  figure  la  leads  to  the  formation  of  positive  Gaussian  curvature, 
whereas  the  pattern  from  figure  1  c  leads  to  negative  Gaussian  curvature.  On  removal  of  the  heat, 
the  film  returns  to  a  largely  flat  state.  It  should  be  noted  that  the  positive  Gaussian  curvature 
sample  exhibits  a  periodic  buckling  around  a  pair  of  oppositely  faced  edges  of  the  film.  This  is 
likely  owing  to  the  relatively  sharp  change  in  director  angle  with  respect  to  the  resolution  of  the 
patterning  technique  near  these  edges  of  the  film.  This  buckling  highlights  the  limitation  on  the 
curvature  that  can  be  achieved  in  nematic  LC  glasses  with  comparatively  small  strains. 

In  order  to  improve  the  quality  of  the  surfaces  that  are  formed  in  stimulated  nematic  glasses, 
a  discoid  subsection  of  the  patterned  films  was  removed  and  exposed  to  stimulus,  as  shown  in 

figure  3.  In  the  case  of  films  encoded  with  either  of  the  identified  patterns,  the  predicted  smooth 
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Figure  4.  (a)  Side-by-side  comparison  of  positive  (left)  and  negative  (right)  Gaussian  curvature  realization  in  actuated 
elastomeric  films.  (, b )  Comparison  of  shape  selection  of  discs  in  the  negative  Gaussian  curvature  case  depending  on  the  size 
of  the  domain.  For  smaller  radii  (3.7  mm  initial  diameter),  a  saddle  shape  is  formed  as  expected  (right).  When  the  radius  of 
the  disc  is  sufficiently  large  (7.1  mm  initial  diameter),  considerably  more  complex  surfaces  with  wavy  edges  are  formed  (left). 
(Online  version  in  colour.) 


curvature  is  realized  upon  stimulation.  Indeed,  shape  selection  of  the  equilibrium  surface  in 
the  positive  curvature  case  seems  to  be  in  remarkable  qualitative  agreement  with  the  predicted 
solution  of  a  spherical  cap.  In  the  negative  curvature  case  too,  the  equilibrium  surface  appears 
largely  consistent  with  the  hyperboloid  saddle  solution  that  is  predicted  to  be  energetically 
favourable  for  a  glassy  film  at  this  scale.  Owing  to  the  tendency  for  buckling  in  areas  of  the  film 
where  the  director  changes  rapidly  with  respect  to  the  resolution  of  the  patterning  process,  it 
should  be  noted  that  the  curvature  cannot  be  increased  by  simply  scaling  the  pattern  to  smaller 
dimensions.  Instead,  higher  strain  materials  are  needed. 

To  facilitate  larger  curvature  realization,  we  prepare  a  comparatively  high  strain  surface- 
alignable  LC  elastomer  with  X  =  0.65  [23].  Using  the  pattern  depicted  in  figure  la,  positive 
Gaussian  curvature  is  encoded  in  the  elastomeric  film.  As  can  be  seen  in  figure  4 a,  the  film 
encoded  for  positive  Gaussian  curvature  forms  part  of  a  sphere  with  a  slightly  elliptical  distortion. 
Figure  4 a  also  shows  a  complexly  buckled  hyperbolic  surface  that  is  formed  when  an  elastomeric 
disc  encoded  to  give  negative  Gaussian  curvature  is  exposed  to  stimulus.  For  smaller  diameter 
films  encoded  with  negative  Gaussian  curvature,  a  classic  saddle  shape  can  be  observed,  as  shown 
in  figure  4b.  The  surface  that  is  formed  by  the  larger  radius  hyperbolic  disc  of  figure  4 b  might 
be  interpreted  as  a  distorted  periodic  Amsler  surface.  Comparatively,  these  deformations  are 
significantly  larger  than  those  observed  for  the  glassy  films  despite  being  more  than  three  times 
as  thick.  Understanding  the  full  spectrum  of  shape  selection  for  films  encoded  with  negative 
Gaussian  curvature  is  an  area  of  ongoing  consideration. 

These  qualitative  explorations  can  be  made  quantitative  if  (i)  the  X  values  can  be  accurately 
known  (results  are  very  sensitive  to  X),  and  (ii)  accurate  imaging  from  several  vantage  points 
is  made,  so  that  shell  shapes  and  dimensions  can  be  accurately  fitted.  This  work,  also  requiring 
further  synthesis,  is  underway. 


3.  Spiral  director  fields 

(a)  Metric  analysis  in  polar  coordinates 

Now,  consider  a  director  field  n  =  n\e1  +  ri2e2  whose  components  are  given  by 


n\(r,6)  =  cos  ^(r,6)  and  ri2(r,0)  =  sin 


(3.1) 


where  (r,0)  are  the  polar  coordinates  on  co,  and  «^(r,0)  =  t/c(rcos0,rsin0)  is  the  alignment  angle 
field  expressed  in  (r,  0)  coordinates.  The  metric  components  with  respect  to  polar  coordinates  are 


7 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


(3.2) 


rspa.royalsocietypublishing.org  Proc.  R.  Soc.  A  472: 20160112 


Downloaded  from  http://rspa.royalsocietypublishing.org/on  March  10,  2017 


Figure  5.  (a)  A  logarithmic  spiral  nematic  pattern  with  a  =  7t/4.  (b)  The  director  field  defined  by  a(r)  =  r/L  and  (c)  the 
resulting  Gaussian  curvature  distribution  K  =  K(r).  (Online  version  in  colour.) 


where  J  is  the  Jacobian  matrix 


J  = 


drXi  deXi 

drX2  3qX2 


(3.3) 


In  terms  of  the  angle  field  W ,  the  metric  components  are 

an  =  A2  +  (A,_2v  -  A2)  sin2(0  - 
ar0  =  aor  =  ^(A_2v  -  X2)  sin  2 (6  -  V) 

and  age  =  r2[X~2v  -  ( X~2v  —  X2)  sin2(0  -  tf')]. 

By  Gauss's  theorem,  the  Gaussian  curvature  is  given  by 


(3.4) 


K  =  ~—(drr?e  -  3 er?r  +  r;errer  -  r;rr?e  +  r, -  r?rn 


-iQ  t^O 


- e  re  > 


arr 


n 1  ee>- 


(3.5) 


For  nematic  angle  fields  of  the  form  ^{r,0)  =  0  +  a(r),  the  metric  components  in  polar 
coordinates  depend  only  on  the  radius  r,  and  the  Gaussian  curvature  is  given  by 


K  = 


(a"  +  si 


sin(2o()  +  2  a2  cos(2o() 


(3.6) 


The  director  field  takes  the  form  n  =  cos {0  +  a)  e1  +  sin(0  +  a)  e2  =  cos  a(r)  er  +  sin  a(r)  ee,  where 
er  and  ee  are  the  unit  radial  and  azimuthal  vectors,  respectively.  Note  that  a  is  the  angle  that  the 
director  makes  with  the  radial  direction.  For  logarithmic  spiral  patterns,  this  angle  is  constant,  so 
that,  by  (3.6),  the  Gaussian  curvature  generated  by  a  director  field  whose  integral  curves  consist  of 
logarithmic  spirals  vanishes  everywhere  except  at  the  point  defect  at  the  origin.  Figure  5a  shows 
a  logarithmic  spiral  pattern  with  a  =  7t/4. 

By  prescribing  a  =  a(r),  we  can  obtain  director  fields  that  generate  radially  symmetric 
Gaussian  curvature  distributions  K  =  K(r).  An  example  with  a(r)  =  r/L  is  shown  in  figure  5b. 
The  Gaussian  curvature,  in  this  case,  is  given  by 


m=- 


[pC05(r)  +  Ssln(r).' 


For  constant  K  e  R,  (3.6)  can  be  viewed  as  an  ODE  in  a  which  is  solved  by 


1  r1  ci 

a(r)  =  ±-  arccos  -C(K)—  +  c\  + 

2  \  2  rl 


(3.7) 


(3.8) 


where  C(K)  =K/(X  2  —  X2v),  and  c\,  C2  are  real  constants  of  integration.  Depending  on  the  choice 

of  c\  and  C2,  these  solutions  define  curvature-inducing  nematic  patterns  on  compact  discs  or 
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Figure  6.  Constant-curvature-inducing  nematic  spiral  patterns  on  discs  of  maximal  radius  td  =  2/-v/|C(K)|.  (a)  The 
director  field  of  the  nematic  spiral  pattern  which  generates  prescribed  constant  curvature  K  >  0.  [b)  The  director  field  of  the 
nematic  spiral  pattern  which  generates  prescribed  constant  curvature  K  <  0.  (c)  Orthogonal  duality  of  constant-curvature- 
inducing  spirals  of  opposite  sign  on  a  disc  of  maximal  radius.  The  spiral  pattern  (red)  whose  integral  curves  turn  in  a  clockwise 
direction  away  from  the  origin  generates  negative  Gaussian  curvature.  (Online  version  in  colour.) 


annular  domains.  If  K  >  0,  C2  =  0  and  c\  <  1,  then  the  solution  (3.8)  defines  a  spiral  pattern  on 
the  compact  disc 


r  < 


2(1 +  d) 

C(K) 


(3.9) 


Note  that  on  the  boundary  of  this  disc,  the  expression  for  a  in  (3.8)  attains  its  maximum  value  of 
a  =  n /2.  For  c\  =  1,  we  obtain  a  spiral  pattern  on  a  disc  of  maximal  radius,  which  encodes  constant 
positive  Gaussian  curvature  K.  This  gives  the  largest  disc  on  which  a  smooth  radially  symmetric 
director  field  can  be  used  to  encode  a  prescribed  constant  Gaussian  curvature.  Attempting  to 
increase  the  radius  beyond  this  point  by  varying  c\  as  a  parameter  leads  to  a  bifurcation  at  which 
the  solution  domain  becomes  annular  instead  of  discoid.  If  K  <  0,  then  it  is  C2  =  0  and  c\  =  —  1 
which  yield  the  spiral  pattern  of  maximal  radius.  These  patterns  are  shown  in  figure  6  a  Jo. 

In  the  light  of  the  result  that  the  orthogonal  dual  of  a  director  field  has  precisely  the  opposite 
Gaussian  curvature  distribution,  we  find  that  the  spirals  in  figure  6  make  perfect  sense.  In 
particular,  the  spiral  in  figure  6b  is  equivalent  to  the  spiral  that  turns  in  the  opposite  direction, 
defined  by  an  a(r)  =  —  \  arccos(— C(K)r2/2  —  1).  This  equivalent  pattern  is  the  orthogonal  dual  of 
the  positive  curvature  spiral  of  figure  6 a  and  is  illustrated  in  figure  6c,  noting  C(K)  has  reversed 
sign. 

Finally,  note  that  for  K  =  0  and  the  choice  of  C2  =  0,  solution  (3.8)  reduces  to  a  logarithmic 
spiral  pattern  defined  by  a  =  constant.  It  should  be  noted,  however,  that  the  resulting  surfaces 
upon  stimulation  will  typically  not  remain  planar.  In  particular,  such  patterns  are  expected  to 
form  cone  and  anti-cone  equilibrium  configurations  upon  stimulation  [4,7,26,27]. 


(b)  Logarithmic  spiral  patterns:  cone/anti-cone  formation 

It  is  shown  in  [4,7]  that  azimuthal  (a  =  tt/2)  and  radial  [a  =  0)  director  fields  generate  cone  and 
anti-cone  surfaces  upon  stimulation,  respectively.  It  is  also  argued  there  that  logarithmic  spiral 
patterns  result  in  cone  and  anti-cone  surfaces  depending  on  whether  the  angle  a  is  greater  than 
or  less  than  some  threshold  angle  depending  on  X  and  v.  Here,  we  will  rederive  this  result  using 
the  nematic  metric  in  polar  coordinates.  Generally,  the  angles  between  curves  on  a  patterned 
nematic  sheet  become  distorted  upon  stimulation.  In  particular,  the  radial  curves  on  the  initial 
sheet  will  not  in  general  be  orthogonal  to  the  images  of  the  azimuthal  curves  upon  stimulation. 

We  note  that  two  curves  Tf  and  F \  with  polar  parametric  representations  (r,0)  =  (ri(£),0i(£))  and 
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(r,0)  =  {j2 (t),  02(f))  are  orthogonal  with  respect  to  the  metric  ds2  =  arr  dr2  +  2 are  dr  d 0  +  aee  d 02,  if 
and  only  if 

ttrrY  1^2  "h  ttrOV\02  +  +  aee  0-^0 2  =  0.  (3.10) 

The  concentric  circles  on  the  undeformed  sheet  can  be  parametrized  such  that  =  0  and  0^  =  1. 
Thus,  the  curve  (r(t),0(t))  on  the  undeformed  sheet  whose  image  becomes  orthogonal  to  the 
images  of  the  concentric  circles  upon  stimulation  satisfies  aerr'  +  aeeO'  —  0.  That  is,  dO/dr  = 
—aor/aQQ.  The  length  of  such  a  curve  emanating  from  the  origin  is  given  by 


h  = 


~r~  aee  dr 


upon  stimulation.  Because  d# /dr  =  —aor/aoo,  (3.11)  reduces  to 

\fhr098  —  i_v  rdr 


h  = 


r  \/  ClrrUOO  ~ 
+JOQ0 


dr  =  X 


o  V^ee 


(3.11) 


(3.12) 


The  perimeter  of  the  closed  curve  that  the  circle  of  radius  r  deforms  into  following  exposure  to 
stimulus  is  given  by 

rln 


h  = 


y/aoo(r)  d6  =  2n-Jaee(r). 


(3.13) 


o 


The  analysis  so  far  is  valid  for  all  angle  fields  of  the  form  &  =  0  +  o'(r),  where  aee  =  2v  — 
(A.-2v  —  A2)  sin2  a]  by  (3.4). 

For  logarithmic  spirals,  a  cone  is  expected  to  form  if  2tt/i  >  I2  and  an  anti-cone  is  expected  to 
form  if  2ttI\  <  I2 .  When  a  cone  is  formed,  the  cone  angle  is 


0  =  arcsin[A.-(1+v)  -  (A-(1+y)  -  A.1+v)  sin2  a]. 


(3.14) 


For  2jtIi  =  I2,  the  sheet  is  expected  to  remain  flat  (0  =  jt/2)  upon  stimulation  and  the  threshold 
angle  ac  at  which  this  holds  satisfies  X~2v  —  (X~2v  —  X2)  sin2  ac  =  A.1-v.  That  is. 


ac  =  arcsin 


I  X~2v  —  X1_v 
X -  X 2 


(3.15) 


which  is  in  agreement  with  [4].  We  also  note  that  we  can  choose  the  angle  a  of  a  logarithmic 
director  field  such  that  the  lengths  of  azimuthal  curves  remain  fixed.  This  is  achieved  by  requiring 
yj aee  (t)  =  r,  which  is  satisfied  from  (3.4)  when 


oiQ  =  arcsin 


\-2v 


-1 


X~2v  -  X 2 


(3.16) 


If  ac  <  oiQ,  then  a  cone  is  formed  and  as  a  result  the  symmetry  of  revolution  is  unbroken.  This 
condition  is  equivalent  to  Xl~v  >  1;  i.e.  v  >  1.  In  particular,  a  logarithmic  spiral  director  field  with 
a  —  oiQ  patterned  across  the  surface  of  a  disc  will  keep  the  circular  boundary  fixed,  whereas  the 
interior  of  the  surface  rises  to  form  a  cone  provided  v  >  1.  The  cone  angle  in  this  case  is  0  = 
arcsin(A.v_1). 
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(c)  Spherical  spindles 

Here  we  extend  the  analysis  presented  in  §3b  to  study  the  formation  of  spherical  caps  and 
spherical  spindles  using  the  positive-curvature-inducing  spiral  patterns  defined  by 

a(r)  =  j  arccos(-  \  C(K)r2  +  c),  (3.17) 

for  a  constant  c.  Spherical  spindles  [28,29]  are  surfaces  of  revolution  of  constant  positive  Gaussian 
curvature  K.  The  top  half  of  a  spherical  spindle  can  be  parametrized  as 

r(s,  t)  =  (yi (s)  cos  t,  yi (s)  sin  t,  y2(s)), 

10 
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Figure  7.  (a)  The  top  half  of  a  spindle  of  constant  Gaussian  curvature  K  >  0  and  parameter  p  e  (0, 1  /VK).  The  spindle 
arises  as  the  surface  of  revolution  obtained  by  revolving  the  profile  curve  y  about  the  x2-axis.  (b)  The  director  field  defined 
by  (3.17)  with  K  >  0  and  c  =  1  —  2/(1  +  A1+y).  The  solid  red  curve  indicates  the  circle  of  radius  ro  whose  length  is 
unchanged  by  the  pattern,  (c)  The  director  field  on  a  disc  of  radius  ro  and  [d)  the  spherical  cap  of  fixed  boundary  that  is  expected 
to  form.  (Online  version  in  colour.) 


where  the  profile  curve  y(s)  =  (/i(s),  /2(s))  is  defined  by 

/i(s)  =  psin(sVK) 

and 


Yz(s)  = 


■ n/2-sVK 


0 


—  —  p2  sin2  r  dr, 
K 


(3.19) 


(3.20) 


for  0  <  p  <  1/VK  and  s  e  [0,7r/2 VK].  Figure  7a  illustrates  one  such  surface  and  its  profile  curve. 
Note  that  for  p  =  1/VK,  the  spindle  reduces  to  a  sphere  of  radius  1/VK. 

The  length  of  each  curve  from  the  apex  that  intersects  the  azimuthal  circles  orthogonally  is 


h  (s)  = 


r=s, 


because  y^2  +  y^  =  1.  Using  (3.12)  we  also  have 


h  =  s(r)  = 


o  y/aee(r) 


By  (3.17),  we  find  that  a$o  =  r2{/ii  —  pL^r2),  where 

/xi  =  \[7~lv  +  X2  +  c(k~lv  —  A.2)]  and 
Hence,  (3.22)  can  be  rearranged  and  exploited  to  get  r(s): 


„  1  <\  2— 2v 

M2  =  4A 


K. 


(3.21) 


(3.22) 


(3.23) 


s(r)  = 


1  l-V 


dr 


0  Tmi 


'  P'2' 


1 1-v 


;  arctan 


1  - 


M2r 


(3.24) 


so  that  r(s)  =  Vmi 7m 2  sin(  \  VKs).  The  length  of  the  image  of  the  circle  of  radius  r  centred  at  the 
origin  in  the  undeformed  sheet  becomes  /2  =  TjiryJ \jl\  —  /x2r2  upon  stimulation.  Thus,  we  find  that 


2~;2(s)  =  r(s)J in  -  7s}2  =  Sin sVK. 


(3.25) 


1  v-1 


By  comparing  the  circumference  /2  with  271/1,  the  spindle  radius  being  (3.19),  we  identify 

/i i/VK  with  p.  Given  p  <  1/y/K,  the  deformed  sheet  is  consistent  with  a  spherical  spindle 
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precisely  when  Xv  1  m  <  1.  That  is,  a  spindle  is  expected  to  form  if 


2X1+V  -  A2(1+u>  -  1 

C  <  - ,  r -  =  1  — 


1  -  A.2<1+y)  1  +  Xl+V  ' 

The  angle  of  the  spindle  at  the  apex  is  <p  =  arcsinf  ]/((()))  =  arcsinOu-i/A.1-1’).  That  is, 

1  +  c  -  (c  -  l)X2+2v 


0  =  arcsin 


2A,1+y 


(3.26) 


(3.27) 


If  c  =  1  —  2/(1  +  A,1+y),  then  a  spherical  cap  of  radius  1/VK  will  form  instead. 

For  any  given  choice  of  the  parameter  c,  the  spiral  in  figure  6 a  keeps  the  circumference  of  an 
azimuthal  circle  of  radius  r  =  tq  fixed  provided  that  it  satisfies  y/aee  (r)  =  r;  see  (3.13).  That  is,  tq 

satisfies  yj  fi\  —  \i2r^  =  1/  so  that 

r°  =  =  /4/l  +  c-A2y(2  +  (c-l)A2).  (3.28) 

V  M2  xVk  v 


In  the  case,  c  —  1  —  2/(1  +  A,1+y),  corresponding  to  a  spherical  cap,  we  obtain 


vV-1 


Vh-D. 


(3.29) 


Note  that  we  require  y  >  1  for  ro  to  be  well-defined  and  positive.  We  have  thus  identified  a  director 
field  on  a  thin  disc  that  will  cause  it  to  form  a  spherical  cap  with  the  same  fixed  circular  boundary 
upon  exposure  to  stimulus,  as  shown  in  figure  7b-d.  The  zenith  angle  cpz  of  the  resulting  spherical 
cap  is 

(pz  =  VXs(ro)  =  2  arctan  \A,1-y  -  1.  (3.30) 
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Note  that  even  for  strains  as  small  as  X  =  0.98  and  an  optothermal  Poisson  ratio  of  v  =  2,  the 
zenith  angle  of  the  resulting  spherical  cap  can  be  as  large  as  approximately  16°.  So  we  see  that 
this  method  allows  for  rather  prominent  spherical  caps  to  form  even  in  the  presence  of  the  small 
strains  observed  in  nematic  glasses.  The  same  methods  can  be  used  to  specify  patterns  that 
generate  spindles  with  unchanged  circular  boundaries. 


(d)  Hyperbolic  cones 

We  can  perform  a  similar  analysis  in  the  case  of  the  negative-curvature-inducing  spiral  patterns. 
The  family  of  surfaces  that  form  the  negative  curvature  analogue  of  spherical  spindles  are 
hyperbolic  cones  [28],  which  are  cone-like  surfaces  of  revolution  of  constant  negative  Gaussian 
curvature  away  from  the  non-smooth  tip  of  the  cone.  A  hyperbolic  cone  of  constant  negative 
Gaussian  curvature  K  <  0  admits  a  parametrization 


r(s,  t)  =  (yi(s)  cos  t,  yi(s)  sin  t,  y2(s)), 
where  the  profile  curve  y(s)  =  (yi(s),  y2(s))  is  defined  by 

yi(s)  =  P  sinh(s/jK|) 

and 

Y2 (s)  = 


fSvW  /  1 

— —  p2  cosh2  r  dr 

\K\ 


for  0  <  p  <  1/vW  and 


„  1  .  u  y/l/\K\  ~  P2 

0  <  s  <  - _ arsmh - . 


vm 
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Figure  8.  (a)  A  hyperbolic  cone  and  its  profile  curve  y.  (b,c)  The  director  field  defined  by  (3.17)  with  K  <  0  and  c  =  — 1. 
The  solid  (red)  circle  indicates  the  circle  of  radius  ro  whose  length  is  unchanged  by  the  stimulation  of  the  pattern.  (Online 
version  in  colour.) 


Note  that  y2  can  also  be  expressed  in  terms  of  the  incomplete  elliptic  integral  of  the  second  kind 
E,  as 


~P2E 


ym  -p1)' 


(3.35) 


where 


E(P\k2)  = 


yjl—k2  sin2  r  dr. 


(3.36) 


Figure  8a  illustrates  a  typical  hyperbolic  cone  and  its  profile  curve  arising  from  a  negative 
curvature  spiral. 

Using  the  spiral  pattern  specified  by  (3.17)  with  prescribed  constant  K  <  0,  the  analysis 
proceeds  as  in  the  case  of  positive  Gaussian  curvature.  The  expression  in  (3.24)  is  modified  to 


s(r)  = 


ll-V 


dr 


o  V/xi  +  |/x2| r2 


1  l-V 


arsinh  (  r. 


I  M2 1 
Ml 


(3.37) 
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where  i±\  and  /x2  are  as  in  (3.23)  with  K  <  0.  This  yields  r(s)  =  VM1/IM2I  sinh( \ y/\K\s),  and  thus 

T/2(s)  =  rJEl  +  lM2|r2  =  sinhsTl^f.  (3.38) 

We  see  that  this  coincides  with  (3.32),  using  y\  as  a  spindle  radius  to  generate  a  circumference,  if 
we  identify  7v-1/xi/VfKT  with  p,  whereby  we  note  that  the  deformed  sheet  is  consistent  with  a 
hyperbolic  cone  precisely  if  A.v-1/x2  <  1;  i.e.  if  c  <  1  —  2/(1  +  A,1+v). 

For  any  given  choice  of  the  parameter  c,  the  spiral  in  figure  6b  keeps  the  circumference  of  an 
azimuthal  circle  of  radius  v  =  tq  fixed  provided  that  it  satisfies  yjoee  (r)  =  r.  That  is,  ro  must  satisfy 

/mi  +  l/i2  |Jq  =  r  so  that 

r“=/W=i^i^’a+(c-1)l2)-a+^  <3'3,) 

Note  that  ro  must  also  satisfy  ro  <  rmax(c),  where 

hnax(c)  =  ^^-^2(1  —  c)(A  2  —  l2v)  (3.40) 

is  the  maximum  radius  on  which  (3.17)  is  well  defined  for  a  prescribed  K<  0.  Figure  8b 

shows  the  negative  curvature  spiral  pattern  defined  by  (3.17)  with  c  =  —  1.  The  solid  (red)  circle 
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Figure  9.  Extension  of  the  positive  curvature  encoding  spiral  pattern  using  solutions  in  the  form  of  (3.8)  defined  on  annular 
domains,  (a)  The  angle  a  that  the  directors  make  with  the  radial  direction  plotted  against  the  radius  r.  An  extension  using  two 
annular  solutions  is  shown.  (, b )  A  visualization  of  the  corresponding  nematic  director  field.  (Online  version  in  colour.) 

superimposed  on  the  director  field  indicates  the  circle  of  radius  tq  whose  length  is  unchanged. 
The  pattern  in  figure  8c  is  thus  expected  to  generate  a  hyperbolic  cone  of  fixed  boundary  upon 
stimulation. 

(e)  Annular  extensions  of  curvature-inducing  spiral  patterns 

Finally,  we  note  that  it  is  possible  to  radially  extend  the  constant  curvature  encoding  spiral 
patterns  of  figure  6  to  cover  larger  domains  by  the  use  of  director  fields  defined  on  an  annulus 
whose  inner  radius  matches  that  of  the  disc.  We  can  do  so  consistently  through  judicious  choices 
of  ci  and  C2  ^  0  in  the  solution  (3.8)  for  a,  to  ensure  that  the  directors  agree  on  the  circular 
boundary  between  the  disc  and  the  first  annulus.  We  can  then  repeat  the  process  and  extend  the 
resulting  pattern  further  by  attaching  another  solution  (3.8)  defined  on  a  second  annulus  whose 
inner  radius  matches  the  outer  radius  of  the  first  annulus  using  a  suitable  new  choice  of  ci  and  C2- 
This  process  can  be  continued  to  produce  director  fields  that  encode  constant  Gaussian  curvature 
on  circular  domains  of  any  desired  size.  The  resulting  Gaussian  curvature  is  well  defined  at  all 
points  on  the  domain  with  the  exception  of  the  boundaries  where  the  different  textures  meet. 
Figure  9  depicts  the  nematic  director  field  of  the  extended  pattern  and  shows  how  the  angle 
function  a  =  a(r)  varies  under  such  an  extension  scheme  for  the  K  >  0  spiral.  The  orthogonal 
dual  of  the  extended  pattern  provides  the  corresponding  result  in  the  case  of  negative  Gaussian 
curvature. 

4.  Discussion 

In  general,  surfaces  morphing  to  Gaussian-curved  forms  offer  rich  possibilities,  most  especially 
for  delivering  work  as  microdevices.  Blocking  the  stroke  of  such  morphing  surfaces  creates  stretch 
rather  than  bend.  Normally,  extensile  spontaneous  strains  of  thin  structures,  when  blocked  lead 
to  Euler  instabilities  and  potentially  strong  processes  become  weaker  bend  actions.  A  Gaussian- 
curved  structure  is  less  likely  to  find  such  soft  alternative  routes.  Impeding  motion  of  the  tip  of 
an  emerging  cone,  and  thus  resisting  the  contraction  of  circumferences  lead  to  large  associated 
forces,  and  these  thus  consequently  do  work,  or  perhaps  pump  an  enclosed  fluid.  One  might 
term  this  action  'strong  actuation'  as  opposed  to  a  weak  actuation  associated  with  bend.  Another 
advantage  is  that  which  is  often  a  short  stroke  (associated  e.g.  with  a  4%  length  change  in  a  glass) 

is  larger  in  a  morphing  shell,  for  instance  in  the  rise  of  a  cone. 
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To  avoid  the  nugatory  effects  of  soft  elasticity  relieving  strong  tensile  forces  during  actuation 
by  director  rotation,  one  must  only  have  such  stresses  along  the  director.  Often  this  can  be 
achieved,  for  instance  when  a  cone  rising  on  heating  is  blocked,  the  circumferential  stresses 
arising  cannot  be  negated  by  director  rotation,  because  stresses  are  already  aligned  along  the 
director.  However,  more  complex  director  distributions  may  be  less  reliable  and  it  would  be  safer 
to  use  elastomers  heated  to  isotropy,  or  glasses  where  there  is  no  director  rotation. 

The  anchoring  of  deforming  shells  to  unresponsive  mounts  surrounding  them  is  of  vital 
importance  for  exploiting  shells  with  tuneable  Gaussian  curvature.  We  have  given  explicit  forms 
for  the  director  distribution  to  achieve  deformations,  of  both  signs  of  curvature,  elastically 
compatible  with  rigid  connections  to  a  mounting  at  the  shell  boundary.  Negative  curvature 
surfaces  might  be  of  practical  advantage  because  they  have  sharp  tips  and  are  thus  more  acutely 
sensed,  for  instance,  in  haptic  devices.  Further,  the  slopes  of  their  profile  curves  decrease  with 
increasing  radius,  in  contrast  to  cones.  Thus,  at  their  fixed  boundaries  where  they  could  be 
mounted,  there  is  less  of  a  difference  in  angle  and  thus  lower  elastic  energy  in  the  attachment 
region.  It  is,  however,  impossible  to  design  the  pattern  so  as  to  ensure  that  the  radius  at  which  the 
boundary  is  fixed  precisely  matches  the  radius  at  which  the  profile  curve  fully  flattens. 

In  summary,  our  results  clearly  indicate  that  Gaussian  curvature  can  be  realized  in  both 
low-strain,  high-modulus  glassy  and  high-strain,  low-modulus  elastomeric  LC  solid  surfaces 
using  appropriate,  smooth,  in-plane  director  fields  patterned  across  initially  flat  films.  Our 
preliminary  experimental  investigations  show  that  patterned  nematic  solids  may  indeed  be 
promising  candidates  for  applications.  In  particular,  the  observed  shape  transformations  were 
in  qualitative  agreement  with  theoretical  predictions  and  the  behaviour  of  the  films  in  response 
to  stimulus  was  found  to  be  robust  and  reproducible.  We  hope  that  our  results  will  encourage 
and  stimulate  further  experimental  research  in  achieving  desired  shape  transitions  in  nematic  LC 
solid  surfaces.  In  particular,  our  discovery  of  anchorable  structures  should  lead  to  further  work 
in  assessing  the  viability  of  morphing  shells  for  specific  applications. 

Data  accessibility.  All  experimental  data  in  this  work  (photographs,  experimental  conditions,  etc.)  are  given  in  the 
text.  All  theoretical /computational  data  are  also  given  in  the  text  and  relate  to  plots  of  mathematical  results 
derived  analytically. 
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